{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# PID step responses" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here are some open loop step responses of PID controllers in different configurations." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## PI" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import control\n", "import numpy\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "s = control.tf([1, 0], 1)\n", "ts = numpy.linspace(0, 5)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def plot_step_response(G):\n", " t, y = control.step_response(G, ts)\n", " # Add some action before time zero so that the initial step is visible\n", " t = numpy.concatenate([[-1, 0], t])\n", " y = numpy.concatenate([[0, 0], y])\n", " plt.plot(t, y)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "K_C = 1" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "tau_I = 1" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "Gc = K_C*(1 + 1/(tau_I*s))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAcFElEQVR4nO3deXhU9f328fcHCIRA2MNOCBh2EkHDJloVN1CoArbVWrUuxfax/dnlJwRFpW6IWmvrUgt1rVofSwCRRYQHN1xQcJmELOz7EnYCIet8nz9IW6sIQ5jJOTO5X9fFRSZzMrkPgTtfTs75HHPOISIi/lXH6wAiInJ8KmoREZ9TUYuI+JyKWkTE51TUIiI+Vy8SL9qqVSuXkpISiZcWEYlJK1as2O2cSzrWcxEp6pSUFJYvXx6JlxYRiUlmtvG7ntOhDxERn1NRi4j4nIpaRMTnVNQiIj6nohYR8bmQitrMmpnZDDPLN7M8MxsS6WAiInJUqKfn/Ql4yzl3pZnVBxIimElERL7mhCtqM2sCfA94FsA5V+ac2x/pYCIi0eSzDXt55r21EXntUA59dAV2Ac+b2Rdm9jcza/TNjcxsnJktN7Plu3btCntQERE/OlRawd1v5PCDZz7m1WWbKC6rCPvnCKWo6wFnAH9xzvUHDgOZ39zIOTfNOZfhnMtISjrmVZAiIjHl3YJCLvnj+/z9k43cMDSFBbedQ0L98F/wHcorbgG2OOeWVT2ewTGKWkSktth3uIz75uUy8/OtpLZuzIyfn8WZnZtH7POdsKidczvMbLOZ9XDOFQAXALkRSyQi4lPOOeZn7+CeOTnsLy7nV8NS+eWwVBrUqxvRzxvqGv1XwCtVZ3ysA26IXCQREf8pPFjCpNk5vJ27k7QOTXnpxkH0bt+kRj53SEXtnPsSyIhwFhER33HO8c/lW7hvXi5lFUEmjujJTWd3oV7dmrteMCJjTkVEYsHmvcVMnJnN0jW7GdilBQ+NSaNrUuMaz6GiFhH5hsqg48WPNvDIwgLq1jHuu6Iv1wxMpk4d8ySPilpE5GtW7yxifFaALzbt5/weSTwwOo32zRp6mklFLSIClFUEeea9tTy5ZA2NGtTl8R/14/J+7THzZhX9dSpqEan1Alv2M35GgPwdRYxMb8fk7/ehVeMGXsf6NxW1iNRaJeWV/HHRKqZ/sI6kxAZMvy6Di3q38TrWt6ioRaRW+mTdHjKzAmzYU8zVAzuROaIXTRvGeR3rmFTUIlKrFJWUM/WtfF7+ZBPJLRJ49eZBnJXayutYx6WiFpFaY0n+Tu6clcPOgyXcdHYXfndx94gMUQo3/ycUETlFew+Xce+bK5n95Ta6t2nM09ecRf/kyA1RCjcVtYjELOccbwa2M3nOSopKyrntgm7cen4q9etF1+1iVdQiEpN2HDg6RGlx3k5O79iUqVcOomfbmhmiFG4qahGJKc45XvtsMw/Oy6M8GGTSZb24YWgX6np0+Xc4qKhFJGZs3HOYzKxsPl63hyFdW/LQ2DQ6t/zWnQOjjopaRKJeZdDx/IfrefTtAuLq1GHKmDSuGtDJF5d/h4OKWkSiWsGOo0OUvtq8nwt6tub+0X1p19TbIUrhpqIWkahUVhHk6XfX8NQ7a0iMj+PPV/dnVHq7mFlFf52KWkSizlebjw5RKthZxOX92nP3yN609NEQpXBTUYtI1DhSVsljiwp4dul6WifG8+z1GVzQy39DlMJNRS0iUeGjtbvJzMpm095ifjwomcwRPWkS788hSuGmohYRXztYUs6U+fn849NNpLRM4B8/G8yQ01p6HatGqahFxLcW5+7kztnZ7CoqZdz3uvKbC7vTsH5dr2PVOBW1iPjOnkOlTH4zlze/2kbPtolMvy6D9I7NvI7lGRW1iPiGc445X21j8pyVHCqt4LcXdefn554WdUOUwi2kojazDUARUAlUOOcyIhlKRGqfbfuPMGl2DkvyC+nXqRkPX5lO9zaJXsfyhZNZUZ/vnNsdsSQiUisFg45/fLaJKfPzqQy6mBiiFG469CEinlm/+zCZWQGWrd/L0NSWTBmdTnLLBK9j+U6oRe2At83MAX91zk375gZmNg4YB5CcnBy+hCIScyoqgzz34Xr+8PYq6terw9SxafwwI3aGKIVbqEU91Dm3zcxaA4vMLN859/7XN6gq72kAGRkZLsw5RSRG5G0/yISsAIEtB7iwVxseGN2XNk3ivY7layEVtXNuW9XvhWY2CxgIvH/8jxIR+Y/SikqeWrKGp99dS9OGcTz54/5clhabQ5TC7YRFbWaNgDrOuaKqty8G7o14MhGJGZ9v2seEGQFWFx5iTP8O3DWyN80b1fc6VtQIZUXdBphV9V2vHvCqc+6tiKYSkZhQXFbBIwsLeOGjDbRrEs/zNwzg/B6tvY4VdU5Y1M65dcDpNZBFRGLI0tW7yZwZYMu+I1w3pDPjh/ekcQOdaFYd+lMTkbA6cKScB+fl8X+Xb6ZLq0a8fssQBnZp4XWsqKaiFpGwWbhyB3fNzmHP4TJ+cd5p3HZBN+Ljat8QpXBTUYvIKdtVVMrkOSuZl72dXu2a8Oz1A0jr2NTrWDFDRS0i1eacY9YXW7l3bi7FpZXcfkkPxn2vK3F1a/cQpXBTUYtItWzdf4Q7Zmbz3qpdnNm5OVPHppPaurHXsWKSilpETkow6Hh52UamLsjHAfeM6s11Q1I0RCmCVNQiErK1uw6RmRXgsw37OKdbKx4cnUanFhqiFGkqahE5oYrKINM/WM8fF68ivl4dHrkynSvP7KjLv2uIilpEjmvltgNMyAqQs/Ugw/u05d4r+tA6UUOUapKKWkSOqaS8kieWrOaZ99bRPKE+f7nmDEaktfM6Vq2kohaRb1mxcS/jZwRYu+swY8/oyF0je9EsQUOUvKKiFpF/O1x6dIjSix9voH3Thrx440DO7Z7kdaxaT0UtIgC8v2oXE2dms+3AEa4b3JnbNUTJN/RVEKnl9heXcf+8PGas2ELXpEb885YhZKRoiJKfqKhFarG3crYzafZK9hWXcev5p/GrYRqi5EcqapFaqLCohHveWMmCnB30ad+EF28cQJ/2GqLkVypqkVrEOceMFVu4f14eR8orGT+8Bz87R0OU/E5FLVJLbNlXzMSZ2XywejcDUprz0Nh0TkvSEKVooKIWiXHBoOOljzfw8MICDLj38j78ZFBn6miIUtRQUYvEsDWFRUzIymbFxn2c2z2JB0b3pWNzDVGKNipqkRhUXhlk2vvr+NPi1SQ0qMtjPzyd0f07aIhSlFJRi8SYnK0HuH1GgLztB7ksvR2TR/UhKbGB17HkFKioRWJESXkljy9ezfQP1tGiUX2e+cmZDO/b1utYEgYqapEY8On6vWRmBVi3+zA/yujEHZf2omlCnNexJExCLmozqwssB7Y650ZGLpKIhOpQaQVTF+Tz90820rF5Q16+aRBnd2vldSwJs5NZUd8G5AFNIpRFRE7COwWF3Dkzm+0HS7hhaAq3X9KDhPr6T3IsCumramYdgcuAB4DfRjSRiBzXvsNl3Dc3l5lfbCW1dWOyfnEWZyQ39zqWRFCo334fB8YDid+1gZmNA8YBJCcnn3oyEfkvzjnmZ+/gnjk57C8u53+GpXLrsFQa1NMQpVh3wqI2s5FAoXNuhZmd913bOeemAdMAMjIyXNgSigiFB0uYNDuHt3N3ktahKS/dOIje7XUUsrYIZUU9FPi+mV0KxANNzOxl59xPIhtNRJxz/HP5Fu6bl0tZRZCJI3py09ldqKchSrXKCYvaOTcRmAhQtaL+X5W0SORt2lPMxFkBPlyzh4FdWjB1bDpdWjXyOpZ4QD8iFvGZyqDjhY828OjCAurWMR4Y3ZerByRriFItdlJF7Zx7F3g3IklEhNU7ixifFeCLTfs5v0cSD4xOo32zhl7HEo9pRS3iA2UVQZ55by1PLllDowZ1+dNV/fj+6e01REkAFbWI5wJb9jN+RoD8HUWMOr09k0f1pmVjDVGS/1BRi3ikpLySPy5axfQP1pGU2IDp12VwUe82XscSH1JRi3jgk3V7yMwKsGFPMVcP7MTES3vRJF5DlOTYVNQiNaiopJyHFuTzyrJNJLdI4NWbB3FWqoYoyfGpqEVqyJL8ndw5K4edB0u4+ewu/O7iHjSsr8u/5cRU1CIRtvdwGfe+uZLZX26je5vGPH3NWfTXECU5CSpqkQhxzvFmYDuT56ykqKSc2y7oxq3np1K/ni7/lpOjohaJgB0Hjg5RWpy3k9M7NmXqlYPo2VZDlKR6VNQiYeSc47XPNvPgvDzKg0HuvLQXN57dhbq6/FtOgYpaJEw27jlMZlY2H6/bw5CuLXlobBqdW2qIkpw6FbXIKaoMOp7/cD2Pvl1AXJ06TBmTxlUDOunybwkbFbXIKSjYcXSI0leb93NBz9bcP7ov7ZpqiJKEl4papBrKKoI89c4ann53DYnxcfz56v6MSm+nVbREhIpa5CR9sWkfE7ICrNp5iCv6tefuUX1o0ai+17EkhqmoRUJUXFbBY2+v4rkP19OmSTzP/TSDYT01REkiT0UtEoKP1uwmc2Y2m/YW8+NByUwc0ZNEDVGSGqKiFjmOA0fKmTI/j9c+20xKywReGzeYwV1beh1LahkVtch3WJS7k0mzs9lVVMot53blNxd2Jz5OQ5Sk5qmoRb5h96FSJs9ZydzAdnq2TWT6dRmkd2zmdSypxVTUIlWcc7zx5TZ+/+ZKDpdW8ruLunPLuadpiJJ4TkUtAmzbf4RJs3NYkl9I/+RmPDw2nW5tEr2OJQKoqKWWCwYdr366iYcW5FMZdNw9sjfXn5WiIUriKypqqbXW7z5MZlaAZev3MjS1JVNGp5PcMsHrWCLfcsKiNrN44H2gQdX2M5xz90Q6mEikVFQGeXbpeh5btIr69erw8Nh0fpDRUZd/i2+FsqIuBYY55w6ZWRyw1MwWOOc+iXA2kbDL3XaQCVkBsrce4OLebbjvir60aRLvdSyR4zphUTvnHHCo6mFc1S8XyVAi4VZaUcmTS9bwl3fX0iwhjqd+fAaXprXVKlqiQkjHqM2sLrACSAWecs4tO8Y244BxAMnJyeHMKHJKVmw8OkRpTeEhxpzRgbsu601zDVGSKBJSUTvnKoF+ZtYMmGVmfZ1zOd/YZhowDSAjI0MrbvFccVkFjy5cxfMfraddk3iev2EA5/do7XUskZN2Umd9OOf2m9m7wHAg5wSbi3hm6erdZM4MsGXfEa4d3JkJI3rSuIFOcpLoFMpZH0lAeVVJNwQuBKZGPJlINRwoLueB+bm8vnwLXVs14vVbhjCwSwuvY4mcklCWGO2AF6uOU9cBXnfOzY1sLJGT91bODu56I4e9h8v4xXmncdsF3TRESWJCKGd9BID+NZBFpFp2FR0dojQvezu92zXh+Z8OoG+Hpl7HEgkbHbSTqOWcY+bnW7l3bi5Hyiu5/ZIejPteV+LqaoiSxBYVtUSlrfuPcMfMbN5btYszOzdn6th0Uls39jqWSESoqCWqBIOOV5Zt5KEF+Tjg99/vw7WDO1NHQ5QkhqmoJWqs3XWIzKwAn23YxzndWvHg6DQ6tdAQJYl9KmrxvYrKINM+WMfji1cTX68Oj1yZzpVnaoiS1B4qavG1ldsOMCErQM7Wgwzv05Z7r+hD60QNUZLaRUUtvlRSXskTS1bzzHvraJ5Qn79ccwYj0tp5HUvEEypq8Z3lG/YyISvA2l2HufLMjky6rBfNEjRESWovFbX4xuHSCh5ZWMCLH2+gfdOGvHTjQL7XPcnrWCKeU1GLL7y3ahd3zMxm24EjXD8khdsv6UEjDVESAVTU4rH9xWXcOzeXmZ9v5bSkRvzzliFkpGiIksjXqajFM/Ozt3P3GznsKy7nl+en8sthqRqiJHIMKmqpcYUHS7j7jZW8tXIHfdo34cUbB9KnvYYoiXwXFbXUGOccM1Zs4b65uZRUBJkwvCc/O6cL9TRESeS4VNRSIzbvLeaOWdl8sHo3A1Na8NDYNLomaYiSSChU1BJRwaDjpY838PDCAgy47/I+XDNIQ5REToaKWiJmTWERE7KyWbFxH+d2T+LBMWl0aNbQ61giUUdFLWFXXhnkr++t5c//bw0JDery2A9PZ3T/DhqiJFJNKmoJq5ytB7h9RoC87Qe5LL0dk0f1ISmxgdexRKKailrCoqS8kscXr2b6B+to2ag+f732TC7p09brWCIxQUUtp+zT9XvJzAqwbvdhfpTRiTsu7UXThDivY4nEDBW1VNuh0gqmLsjn759spFOLhrxy8yCGprbyOpZIzFFRS7W8U1DInTOz2X6whBuHduF/L+lOQn39dRKJBP3LkpOy93AZ983NZdYXW0lt3ZisX5zFGcnNvY4lEtNOWNRm1gl4CWgLBIFpzrk/RTqY+ItzjnnZ27nnjZUcOFLO/wxL5dZhqTSopyFKIpEWyoq6Avidc+5zM0sEVpjZIudcboSziU/sPFjCpNk5LMrdSVqHprx88yB6tWvidSyRWuOERe2c2w5sr3q7yMzygA6AijrGOed4fflm7p+XR1lFkIkjenLT2RqiJFLTTuoYtZmlAP2BZcd4bhwwDiA5OTkM0cRLm/YUkzkzwEdr9zCwSwumjk2nS6tGXscSqZVCLmozawxkAb92zh385vPOuWnANICMjAwXtoRSoyqDjhc+2sCjCwuoW8d4YHRfrh6QrCFKIh4KqajNLI6jJf2Kc25mZCOJV1bvLGJ8VoAvNu1nWM/WPDC6L+2aaoiSiNdCOevDgGeBPOfcY5GPJDWtrCLIM++t5Yklq2ncoB6P/6gfl/drryFKIj4Ryop6KHAtkG1mX1a97w7n3PzIxZKa8tXm/UzICpC/o4hRp7dn8qjetGysIUoifhLKWR9LAS2tYsyRskoeX7yK6R+sIymxAdOvy+Ci3m28jiUix6ArE2uhT9btITMrwIY9xVw9sBMTL+1Fk3gNURLxKxV1LXKwpJyHFuTz6rJNJLdI4NWbB3GWhiiJ+J6KupZYkr+TO2bmUFhUws/O6cJvL+pBw/q6/FskGqioY9yeQ6XcOzeXN77cRo82iTxz7Zn069TM61gichJU1DHKOcebge1MnrOSopJyfn1hN/7PeanUr6fLv0WijYo6Bm0/cIS7ZuewOK+Q0zs14+Gx6fRom+h1LBGpJhV1DAkGHa99tpkp8/MoDwaZdFkvbhjahbq6/FskqqmoY8SG3YfJnBngk3V7Gdy1BQ+NSSdFQ5REYoKKOspVBh3PLV3PHxYVEFenDlPGpHHVgE66/Fskhqioo1j+joNMmBHgqy0HuLBXa+6/Io22TeO9jiUiYaaijkKlFZU89c5ann5nDU0bxvHE1f0Zmd5Oq2iRGKWijjJfbNrHhKwAq3Ye4op+7bl7VB9aNKrvdSwRiSAVdZQoLqvgD2+v4rkP19O2STzP/TSDYT01REmkNlBRR4GP1uwmc2Y2m/YWc82gZDJH9CRRQ5REag0VtY8dOFLOlPl5vPbZZlJaJvDauMEM7trS61giUsNU1D719sodTJqdw+5Dpdxybld+c2F34uM0REmkNlJR+8zuQ6VMnrOSuYHt9GybyN+uzyC9o4YoidRmKmqfcM4x+8ut/P7NXIpLK/ntRd35+bmnaYiSiKio/WDb/iPcOSubdwp20T/56BClbm00RElEjlJReygYdLzy6SamLsinMui4e2Rvrj8rRUOUROS/qKg9sn73YSZkBfh0/V7OTm3FlDFpdGqR4HUsEfEhFXUNq6gM8rel6/njolXUr1eHh8em84OMjrr8W0S+k4q6BuVtP8j4GQGytx7g4t5tuO+KvrRpoiFKInJ8KuoaUFpRyZNL1vCXd9fSLCGOp358BpemtdUqWkRCcsKiNrPngJFAoXOub+QjxZYVG48OUVpTeIgxZ3Tgrst601xDlETkJISyon4BeBJ4KbJRYktxWQWPLCzghY820L5pQ164YQDn9WjtdSwRiUInLGrn3PtmlhL5KLHjYEk5o55YysY9xVw3pDPjh/ekcQMdZRKR6glbe5jZOGAcQHJycrheNiqt2lHExj3FPDw2nR8O6OR1HBGJcmG7Ptk5N805l+Gcy0hKSgrXy0a1ds10RoeInDoNkhAR8TkVdRg453DOeR1DRGLUCYvazP4BfAz0MLMtZnZT5GNFH50TLSKREspZH1fXRJBoppIWkUjSoQ8REZ9TUYuI+JyKWkTE51TUIiI+p6IWEfE5FbWIiM+pqEVEfE5FLSLicypqERGfU1GLiPicilpExOdU1CIiPqeiFhHxORW1iIjPqahFRHxORS0i4nMqahERn1NRi4j4nIpaRMTnVNQiIj6nohYR8TkVtYiIz6moRUR8TkUtIuJzIRW1mQ03swIzW2NmmZEOJSIi/3HCojazusBTwAigN3C1mfWOdDARETmqXgjbDATWOOfWAZjZa8DlQG64w4x6Yikl5ZXhftkadyQG9kFE/COUou4AbP7a4y3AoG9uZGbjgHEAycnJ1QpzWlIjyiqD1fpYvxnctSXpHZp5HUNEYkAoRW3HeJ/71jucmwZMA8jIyPjW86F4/Kr+1fkwEZGYFsoPE7cAnb72uCOwLTJxRETkm0Ip6s+AbmbWxczqA1cBcyIbS0RE/uWEhz6ccxVm9ktgIVAXeM45tzLiyUREBAjtGDXOufnA/AhnERGRY9CViSIiPqeiFhHxORW1iIjPqahFRHzOnKvWtSnHf1GzXcDGan54K2B3GON4KVb2JVb2A7QvfhQr+wGnti+dnXNJx3oiIkV9KsxsuXMuw+sc4RAr+xIr+wHaFz+Klf2AyO2LDn2IiPicilpExOf8WNTTvA4QRrGyL7GyH6B98aNY2Q+I0L747hi1iIj8Nz+uqEVE5GtU1CIiPufLojazH5jZSjMLmlnUnbYTKzcDNrPnzKzQzHK8znKqzKyTmb1jZnlVf7du8zpTdZhZvJl9amZfVe3H773OdKrMrK6ZfWFmc73OcirMbIOZZZvZl2a2PJyv7cuiBnKAMcD7Xgc5WTF2M+AXgOFehwiTCuB3zrlewGDg1ij9upQCw5xzpwP9gOFmNtjjTKfqNiDP6xBhcr5zrl+4z6X2ZVE75/KccwVe56imf98M2DlXBvzrZsBRxzn3PrDX6xzh4Jzb7pz7vOrtIo4WQwdvU508d9ShqodxVb+i9owAM+sIXAb8zessfubLoo5yx7oZcNQVQiwzsxSgP7DM2yTVU3Wo4EugEFjknIvK/ajyODAeiIW7WjvgbTNbUXWz77AJ6cYBkWBmi4G2x3jqTufcGzWdJ4xCuhmweMPMGgNZwK+dcwe9zlMdzrlKoJ+ZNQNmmVlf51zU/RzBzEYChc65FWZ2ntd5wmCoc26bmbUGFplZftX/Sk+ZZ0XtnLvQq88dYboZsE+ZWRxHS/oV59xMr/OcKufcfjN7l6M/R4i6ogaGAt83s0uBeKCJmb3snPuJx7mqxTm3rer3QjObxdHDoGEpah36CD/dDNiHzMyAZ4E859xjXuepLjNLqlpJY2YNgQuBfG9TVY9zbqJzrqNzLoWj/06WRGtJm1kjM0v819vAxYTxm6cvi9rMRpvZFmAIMM/MFnqdKVTOuQrgXzcDzgNej9abAZvZP4CPgR5mtsXMbvI60ykYClwLDKs6ferLqpVctGkHvGNmAY4uChY556L6tLYY0QZYamZfAZ8C85xzb4XrxXUJuYiIz/lyRS0iIv+hohYR8TkVtYiIz6moRUR8TkUtIuJzKmoREZ9TUYuI+Nz/B6PnxdR3ErIiAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot_step_response(Gc)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## PID" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Because the ideal PID is unrealisable, we can't plot the response of the ideal PID, but we can do it for the realisable ISA PID.\n", "\n", "$$G_c = K_C\\left( 1 + \\frac{1}{\\tau_I s} + \\frac{\\tau_D s}{\\alpha \\tau_D s + 1} \\right)$$" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "alpha = 0.1\n", "tau_D = 1" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "Gc = K_C*(1 + 1/(tau_I*s) + 1*s/(alpha*tau_D*s + 1))" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAbFklEQVR4nO3de3CU15nn8e+jbt1aCJBA6IZk7mCMjXBkx/EVxzdsY0RSmayzlalsyjVUxcmMk3iSOJOZzUztTu1Uze5MUrX7D7VJZraSstebi4QvYxvf4niS2AarxR18B7oFEnehu9Rn/+gWMA4eQ7+NTr+t36fKJanVUj+vkX4cznvOecw5h4iIhE+R7wJERCQ7CnARkZBSgIuIhJQCXEQkpBTgIiIhFZ3MF5s9e7abN2/eZL6kiEjobd269YhzrubDj09qgM+bN48tW7ZM5kuKiISemX1wvsc1hSIiElIKcBGRkFKAi4iElAJcRCSkFOAiIiGlABcRCSkFuIhISCnAJ9kLuw/TfXLQdxkiUgAU4JPswZ+9yU9/f941+SIiF0UBPsnGU46+oTHfZYhIAVCAezAwMu67BBEpAApwDwZGNAIXkeAU4B70D2sELiLBKcA90AhcRHJBAe6B5sBFJBcU4B4owEUkFxTgHvQPawpFRIJTgHugEbiI5IIC3IOBkTGcc77LEJGQU4B7kHIwPJbyXYaIhJwC3BPNg4tIUApwTzQPLiJBKcA9UYCLSFAfG+Bm9mMz6zGzHec8Vm1mm83srczbqktbZuHp125MEQnoQkbg/wSs+dBjjwAvOOcWAy9kPpaLMKDzUEQkoI8NcOfcK8CxDz3cBvxz5v1/BtbnuK6Cp/NQRCSobOfAa51z3QCZt3M+6olmtsHMtpjZlt7e3ixfrvBoDlxEgrrkNzGdcxudc63OudaamppL/XKhoTlwEQkq2wA/bGb1AJm3PbkraWoY1AhcRALKNsA3AV/KvP8loCM35UwdauogIkFdyDLCR4HfAUvN7KCZPQD8HXCHmb0F3JH5WC6CbmKKSFDRj3uCc+4LH/Gp23Jcy5SiOXARCUo7MT3RKhQRCUoB7ok28ohIUApwTzSFIiJBKcA90TJCEQlKAe5JvwJcRAJSgHuiZYQiEpQC3BNt5BGRoBTgngxqBC4iASnAPRkYHSeVUmd6EcmeAtyDaJHhHAyNaRpFRLKnAPcgVhIBtBtTRIJRgHsQK0kfQaPdmCIShALcg1hpegSu3ZgiEoQC3IOKiRG4AlxEAlCAe1CuOXARyQEFuAcVmQDXZh4RCUIB7kGsVFMoIhKcAtyDMyNwTaGISAAKcA8mlhFqO72IBKEA96CsWHPgIhKcAtyDSBGUF0c0By4igSjAPYmVRLSMUEQCUYB7EitVgItIMApwTypKovQPawpFRLKnAPdEUygiEpQC3JNYSVQ3MUUkEAW4JxqBi0hQgQLczL5hZjvNbIeZPWpmZbkqrNBVlEZ1nKyIBJJ1gJtZI/BnQKtzbgUQAe7PVWGFLlYSYVAjcBEJIOgUShQoN7MoEAOSwUuaGmIlEe3EFJFAsg5w51wC+O/AfqAbOOmce+7DzzOzDWa2xcy29Pb2Zl9pgYmVRBkcHWdcnelFJEtBplCqgDZgPtAAVJjZFz/8POfcRudcq3OutaamJvtKC0xFpq3a4KhG4SKSnSBTKLcD7znnep1zo8AvgetzU1bhK1dbNREJKEiA7weuM7OYmRlwG7A7N2UVvokzwdWZXkSyFWQO/DXg58CbwPbM99qYo7oK3sSZ4FpKKCLZigb5Yufc94Hv56iWKeXMHLiWEopIlrQT05OY2qqJSEAKcE8mplAGdCKhiGRJAe5JxZk5cI3ARSQ7CnBPyksm5sA1AheR7CjAPZm4iakRuIhkSwHuSVk0gpnmwEUkewpwT4qKjFixzgQXkewpwD0qL4lqCkVEsqYA96iiNKKzUEQkawpwj2IlUZ0JLiJZU4B7FCuJMDiqEbiIZEcB7pG68ohIEApwjypKopoDF5GsKcA9ipVqBC4i2VOAe5SeA1eAi0h2FOAeVZRE6ddOTBHJkgLco1hJlOGxlDrTi0hWFOAeTRxopRuZIpINBbhHE0fK6jwUEcmGAtyjM00dNA8uIllQgHsU0whcRAJQgHt0pi+mAlxEsqAA9yh2piuPplBE5OIpwD2qONOZXiNwEbl4CnCPzs6BawQuIhdPAe6RbmKKSBAKcI8qSjPLCDUCF5EsBApwM5tpZj83sz1mttvMPpWrwqaC0mgRRaY5cBHJTjTg1/8QeMY59zkzKwFiOahpyjAzYiVRTaGISFayHoGb2XTgZuBHAM65EefciVwVNlXEStTYWKRQnRwc5fE3DvDln7x+SX7Pg4zAFwC9wE/MbCWwFXjIOdd/7pPMbAOwAaC5uTnAyxWmitIo/RqBixSModFxXt7bQ3tnkhf39jAylmLerBgHjg2ytK4yp68VJMCjwNXAnzrnXjOzHwKPAH917pOccxuBjQCtra06N/VDYiURBjUCFwm18ZTjtXeP0hFP8vSObvqGxpg9rYT/eG0z61c1snLuDMws568bJMAPAgedc69lPv456QCXi6DGxiLh5JxjZ/IU7Z0JntiW5PCpYSpKIty1oo71LY1cv3AW0cilXeiXdYA75w6Z2QEzW+qc2wvcBuzKXWlTQ6wkyomBEd9liMgF+uBoPx3xJB3xBO/09lMcMW5ZMoe/WtvAbctqzxwTPRmCrkL5U+BnmRUo7wJfDl7S1FJRGiFxQiNwkXx25PQwT23rpj2eoHN/eq3GtfOreeDGBdxzZR0zYyVe6goU4M65ONCao1qmpPLiKIO6iSmSd/qHx3hu1yHaO5O8+vYRxlOOZXWVPHL3Mu5b2UDjzHLfJQYegUtAFaUR7cQUyRMjYyl+81Yv7fEkm3cdYmg0RePMcjbcvIC2lgaW1U33XeK/oQD3LFYS1U5MEY9SKceWD47TEU/w9PZujg+MUhUr5nOfmEtbSyOfaK6iqCj3K0hyQQHuWUVJhJHxFKPjKYov8R1rETlrz6FTdMSTbIonSZwYpKy4iDuX19HW0sBNi2soieb/76MC3LNzGxvPKM//HxiRMEucGGRTZgXJnkN9RIqMmxbP5s/vWsKdy+vOHDAXFuGqtgBN/MAMjIwxo7zYczUihefEwAhPbe+mozPJ6+8fA+Dq5pn8zboruPeqemZPK/VcYfYU4J5NnAmuzTwiuTM4Ms7zuw/TEU/w6329jI47FtZU8PAdS2hraaR5VmGcu6cA92yirZqWEooEMzae4l/fOUpHZ4Jndx6if2Sc2uml/Kfr59HW0sgVDdMvyXZ2nxTgnp0ZgWspochFc84RP3CCjniSJ7clOXJ6hMqyKGuvaqBtVQOfnD+LSJ6uIMkFBbhnsXPmwEXkwrzbe5r2eJJN8QTvHx2gJFrEbcvm0NbSyK3LaiiNTt52dp8U4J5VaA5c5IL0nBpiU1eSjniS7YmTmMH1C2fx4OpF3LWibkouAlCAezaxjFBz4CJ/6NTQKM/sOERHPMHv3jlKysGVjTP4y3sv576VDdROL/NdolcKcM8mbmJqDlwkbXhsnJf29LKpK8Hzu9MNES6bFeNrty5iXUsji+ZM811i3lCAexYrPbuRR2SqSqUcv3/vKJviSZ7e3s2pcxoitLU00NI0s+BWkOSCAtyzkkgR0SLTTUyZciYaImzqSm9nP3RqKN0Q4Yo61rU0cOOi2Ze8IULYKcA9MzPK1ZVHppD9RwfoiCfo6Eryds9pokXG6qU1fO/ey7n98sltiBB2CvA8UFES1QhcCtpEQ4SOeII3JxoizKvmbz+zgntW1FNV4achQtgpwPNArDSizvRScPqHx9i86zDt8QS/eetsQ4TvrFnGupb8aIgQdgrwPFBZVszJgVHfZYgENjqe4pV9vXTEk2zedZjB0fEzDRHWtzSytK7Sd4kFRQGeB+ZWlbMzcdJ3GSJZSaUcW/enGyI8tS3dEGFmrJjPXt1IW0sjrZflb0OEsFOA54GmqhjP7TzEeMoV9LkNUlj2He6jvTNBxzkNEW6/vJb1LY3cvCQcDRHCTgGeB5qrY4yOOw6dGtK8oOS15IlBNnUlae882xDhxkWzefjOJdx5RR3TQtYQIez0fzsPNFWnQ/vAsQEFuOSdEwMjPL09vZ399feP4Rysap7JX9+3nHuvaqCmMrwNEcJOAZ4HmqvTh8vvPzbAdQtmea5GBIZG0w0R2juT/HpfD6PjjgU1FXzj9iW0tTRw2awK3yUKCvC80DCznCKDg8cGfJciU9jYeIrfvnOU9niCZ3ekGyLMqSzlS5+ax/pVhdkQIewU4HmgOFJE/Yxy9ivAZZI55+g6eJL2zgRPbuvmyOlhKkuj3HtVPW0tjVy3oLAbIoSdAjxPNFWXc+D4oO8yZIr4g4YIkSJuXVbD+pZGbl02h7JibWcPAwV4nmiujvHS3l7fZUgB6zk1xBOZ7ezbDqYbIlw3fxZfWb2QNSvqp2RDhLALHOBmFgG2AAnn3NrgJU1NTVUxevuGGRwZ12E+kjN9ZxoiJPntO0dIObiiYTrfuyfdEKFuxtRuiBB2uRiBPwTsBqbn4HtNWc2z0itRDh4fYHGtthtL9obHxnl5by8d8bMNEZqrY3z11kW0tTSwaI5+vgpFoAA3s7nAvcDfAt/MSUVT1NyqdIAfUIBLFlIpx2vvHaMjnjjTEGFWRQlfuKaJtlWNrFJDhIIUdAT+A+DbwEcmjpltADYANDc3B3y5wnV2M49uZMqFcc6xq/sUHfEkT3Ql6T45RCzTEKFNDRGmhKwD3MzWAj3Oua1mtvqjnuec2whsBGhtbXXZvl6hq5lWSllxkZYSysc6cCzTECGe5K1MQ4RbltTwyN3LuGN5LbESrU2YKoL8Sd8ArDOze4AyYLqZ/dQ598XclDa1mBlNVTEOKMDlPI6eHuap7d10xJNs/eA4ANfMq+K/rl/BPVfWU62GCFNS1gHunPsu8F2AzAj8zxXewTRXxzQClzMGRjINETrTDRHGUo6ltZV8e81S1q1sOHPfRKYu/VsrjzRVx3jtvWM453TDaYoaHU/x6ltHaI8neG5nuiFCw4wyHrhpPutbGrm8Xou95KycBLhz7mXg5Vx8r6msqTrG6eExjg+M6p/EU4hzjjf3H6e9M8lT27s51j/CjPJi1q9qpK2lgWvnVashgpyXRuB5pKnq7LGyCvDC99bhPtozNyMPHh+kNFrEHctraWtp5BY1RJALoADPIxObefYfG2Bl00zP1cil0H1ykCe6krR3JtnVfYoigxsWzeYbty/hrhVqiCAXRz8teaTpnM08UjhODozy9I70GSTpexywsmkm379vOWvVEEECUIDnkYrSKLMqSrSUsAAMjY7zwu4eOuIJXt7by8h4igWzK/j6bUtY19LA/NlqiCDBKcDzzNzqmHZjhtR4yvHbd47QEU/yzI5DnB4eY05lKV+87jLWr2rgysYZWl0kOaUAzzPN1TG6DpzwXYZcIOcc2w6eTG9n35akty/dEOHuFXWsX6WGCHJpKcDzTFNVOf+yvZux8ZTOschj7x/ppz2eYFM8ybtH+tUQQbxQgOeZ5uoYYylH98khmqq10y6f9PQN8WRX+mZk1zkNETbcvIC7V9QzI6aGCDK5FOB5ZiK0DxwfUIDngb6hUZ7deZiOeIJ/fftsQ4S/uGcZ961soH5Gue8SZQpTgOeZiaWEB48NwkLPxUxRI2MpXt7bQ0c8yfO7DzM8lqKpupwHVy9i/So1RJD8oQDPM/Uzy4gUmQ61mmSplOON94/RHk/y9PZuTg6mjzP4D9c00dbSyNXNaogg+UcBnmeKI0XUzyjTZp5J4Jxjd3cfHV0JnognSZ4corw4wl1X1NK2qpEbF82mWDeSJY8pwPOQjpW9tA4cG2BTV5KOeIJ9h9MNEW5eUsN31BBBQkY/qXmoqSrGC3t6fJdRUI71j6QbInQm2JJpiNB6WRX/Zf0K7lVDBAkpBXgeap4V48jpYQZGxjQaDGCiIUJHPMkr+3oZSzkWz5nGt+5KN0TQKh8JO6VDHpqbOVb24PFBlqhD/UUZHU/x6ttH6OhM8NyuwwyMjFM/o4wHbpxPW0sjl9dX6makFAwFeB5qzowM9x8dUIBfgHRDhBN0xBM8ta2bo/0jTC+L0tbSQFtLoxoiSMFSgOehczfzyEd7u6eP9s4kHV0JDhxLN0S4fXktbSsbuGVpDaVRbWeXwqYAz0OzKkqIlUS0EuU8Dp0cYlNX4g8aIjx02xLuuqKWyjJtZ5epQwGeh8yMpiodKzvh5OAoz+zopr0zye/fO3qmIcJ/XructSvrmVNZ5rtEES8U4HmqqTrG+0f7fZfhzdDoOC/uSTdEeGlPuiHC/NkVPHTbYtpaGtUQQQQFeN66Zl4Vz+8+zNs9fVPm7I3xlON37xylI57gmR2H6BseoybTEKGtpYGr5qohgsi5FOB56rNXz+Xvn93LY68f4C/XLvddziXjnGNH4hTt8QRPdCXp6RtmWmmUNSvqWN/SyKcWqiGCyEdRgOepmspS7lheyy/ePMi31iwtuBUV7x/ppyOeXkHybm8/xRHj1qVzWL+qkU+rIYLIBVGA57H7r23mX3Yc4rmdh7lvZYPvcgLr7RvmyW1J2uPJM23jPjm/mj+5aQH3qCGCyEVTgOexmxbNpnFmOY+9sT+0AX56eIxndxyi/ZyGCMvrp/Pdu9MNERpmqiGCSLYU4HmsqMi4/5om/sfmfXxwtJ/LZoVj5cXIWIpf7+ulPZ7g+V3phghzq8r5yuqFrG9pZLF2l4rkRNYBbmZNwP8B6oAUsNE598NcFSZpf9TaxD8+v4/H3jjAd9Ys813OR/qohgifb21i/aoGrm6u0goSkRwLMgIfAx52zr1pZpXAVjPb7JzblaPaBKibUcanl83h/205yDfvWJJXDQbONESIJ9jUlaT75BCxkgh3LK9lfUsjNy5WQwSRSynrAHfOdQPdmff7zGw30AgowHPs/muaeX73Fl7Y3cOaFXW+y2H/0QE2dSXoiCd5qyfdEOGWJTU8ooYIIpMqJ79pZjYPWAW8dp7PbQA2ADQ3N+fi5aac1UtrqJ1eymNv7PcW4L19wzy1LUlHV5LO/ekVJNfOq1ZDBBGPAge4mU0DfgF83Tl36sOfd85tBDYCtLa2uqCvNxVFI0V8vrWJ//nS2xw8PsDcqslpRHBycJRndx5iUzzJb99JryBZVlfJI5kVJI1aQSLiVaAAN7Ni0uH9M+fcL3NTkpzPRIA/npkLv1T6h8d4cU8PT3QleXlv+gyS5uoYD65exH0rG1hapxUkIvkiyCoUA34E7HbO/UPuSpLzaaqOcdPiGv7vG/v5fOvcnI7CB0fGeWlvD09uS/Linh6GRlPMyZxBsq6lgZU6g0QkLwUZgd8A/DGw3czimcf+wjn3dPCy5Hy+dusivvyT11nzg9/w/fuW87lPzM06WE8MjPDS3h427zrMS3t6GRwdZ/a0Uj7f2sS9V9bTOq9aZ5CI5Lkgq1BeBfQbPomunV/NM1+/mYcf7+JbP9/Gc7sO898+eyWzp5V+7Nc653int5+XM6G95YPjjKccNZWlfObqRtZeVc8n5+vgKJEw0XqvkGmqjvHohuv48avv8ffP7eXOf3yFv153Bcvrp1MSKSIaMaKRdAjvPdTHmx+coPPAceIHTnBiYBSApbWVfOWWhdy+vJarGmeoX6RISCnAQyhSZPzJzQu4ZWkN33w8zp892vmRzzWDxXOmcdfyOlY1z+T6hbNpnjU5q1hE5NJSgIfYktpKfvXgDbyyr5fTw2OMjjvGxlOMphyplGNhzTSuaprBdPWJFClICvCQK44Ucdvltb7LEBEPdFCFiEhIKcBFREJKAS4iElIKcBGRkFKAi4iElAJcRCSkFOAiIiGlABcRCSkFuIhISCnARURCSgEuIhJSCnARkZBSgIuIhJQCXEQkpBTgIiIhpQAXEQkpBbiISEgpwEVEQkoBLiISUgpwEZGQUoCLiISUAlxEJKQU4CIiIRUowM1sjZntNbO3zeyRXBUlIiIfL+sAN7MI8L+Au4HlwBfMbHmuChMRkX9fNMDXXgu87Zx7F8DMHgPagF25KOxc3/vVdl5/71iuv60XYynnuwQRKRBBArwROHDOxweBT374SWa2AdgA0NzcnNULNcwsZ3HttKy+Nt8sqavk7hX1vssQkQIQJMDtPI/9wfDSObcR2AjQ2tqa1fDzq7cuyubLREQKWpCbmAeBpnM+ngskg5UjIiIXKkiAvwEsNrP5ZlYC3A9syk1ZIiLycbKeQnHOjZnZ14BngQjwY+fczpxVJiIi/64gc+A4554Gns5RLSIichG0E1NEJKQU4CIiIaUAFxEJKQW4iEhImXOTt7XbzHqBD7L88tnAkRyW41OhXEuhXAfoWvJVoVxL0Ou4zDlX8+EHJzXAgzCzLc65Vt915EKhXEuhXAfoWvJVoVzLpboOTaGIiISUAlxEJKTCFOAbfReQQ4VyLYVyHaBryVeFci2X5DpCMwcuIiL/VphG4CIicg4FuIhISIUqwM3sj8xsp5mlzCx0S4sKpQm0mf3YzHrMbIfvWoIysyYze8nMdmd+th7yXVO2zKzMzF43s67MtfyN75qCMLOImXWa2ZO+awnCzN43s+1mFjezLbn83qEKcGAH8FngFd+FXKwCawL9T8Aa30XkyBjwsHPucuA64Ksh/nMZBj7tnFsJtABrzOw6zzUF8RCw23cROXKrc64l12vBQxXgzrndzrm9vuvI0pkm0M65EWCiCXToOOdeAQqiy7Rzrts592bm/T7SgdHot6rsuLTTmQ+LM/+FcpWCmc0F7gX+t+9a8lmoAjzkztcEOpRBUajMbB6wCnjNbyXZy0w7xIEeYLNzLqzX8gPg20DKdyE54IDnzGxrpsl7zgRq6HApmNnzQN15PvU951zHZNeTQxfUBFr8MLNpwC+ArzvnTvmuJ1vOuXGgxcxmAr8ysxXOuVDdqzCztUCPc26rma32XU8O3OCcS5rZHGCzme3J/Cs2sLwLcOfc7b5ruETUBDpPmVkx6fD+mXPul77ryQXn3Akze5n0vYpQBThwA7DOzO4ByoDpZvZT59wXPdeVFedcMvO2x8x+RXo6NScBrimUyaMm0HnIzAz4EbDbOfcPvusJwsxqMiNvzKwcuB3Y47eqi+ec+65zbq5zbh7p35MXwxreZlZhZpUT7wN3ksO/UEMV4Gb2GTM7CHwKeMrMnvVd04Vyzo0BE02gdwOPh7UJtJk9CvwOWGpmB83sAd81BXAD8MfApzPLvOKZkV8Y1QMvmdk20gOGzc65UC/BKwC1wKtm1gW8DjzlnHsmV99cW+lFREIqVCNwERE5SwEuIhJSCnARkZBSgIuIhJQCXEQkpBTgIiIhpQAXEQmp/w+vskAF6d63AQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot_step_response(Gc)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## PD" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "Gc = K_C*(1 + 1*s/(alpha*tau_D*s + 1))" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAARjUlEQVR4nO3dfWxd9X3H8c/H99qOrwPEaQwkMWnalSJaOh7mtnRIXQttRR8E27RNVKLrSrtoKlvpVKlqVWnV/pu0qeo0Tauilj4MRlUVqtKHQRkPRWgUcAKUUMPKaIGQQMxDgMTGju3v/vB1ErJAfM/vJL9zrt8vKbrP93wPhM/98Tvfc36OCAEA6qcndwEAgGIIcACoKQIcAGqKAAeAmiLAAaCmmsdyY2vWrImNGzcey00CQO1t2bLlmYgYPvT5YxrgGzdu1NjY2LHcJADUnu3HDvc8UygAUFMEOADUFAEOADVFgANATRHgAFBTBDgA1BQBDgA1RYAfYzePP62dL0zlLgNAFyDAj7FPX71VV/3isD35ANARAvwYm5sPvfTybO4yAHQBAjyDyZm53CUA6AIEeAaTM4zAAaQjwDNgBA6gDAR4BpPTBDiAdAR4BpP7mEIBkI4Az4AROIAyEOAZMAcOoAwEeAZ76UIBUAICPIOpmTlFRO4yANQcAZ7B7HxoZm4+dxkAao4Az2SKeXAAiQjwTPYS4AASEeCZTE5zIBNAmiMGuO0rbe+yve2g51bbvsn2r9u3Q0e3zO5DKyGAVEsZgX9L0oWHPPcFSTdHxKmSbm4/RgdoJQSQ6ogBHhG3S3rukKcvlvTt9v1vS/rDkuvqehzEBJCq6Bz4SRGxU5Latye+2httb7I9ZntsYmKi4Oa6DwcxAaQ66gcxI2JzRIxGxOjw8PDR3lxtTDGFAiBR0QB/2vZaSWrf7iqvpOVhLxe0ApCoaIBfL+nj7fsfl/TDcspZPqb2EeAA0iyljfAaSXdKOs32dtuflPQPkt5v+9eS3t9+jA7spQ8cQKLmkd4QER99lZcuKLmWZYU+cACpOBMzExY2BpCKAM+EETiAVAR4JgQ4gFQEeCZMoQBIRYBnwggcQCoCPBMCHEAqAjwTrgcOIBUBngkXswKQigDPhMvJAkhFgGcyMzevfaxMDyABAZ7Bit6Ff+wcyASQggDPYLBv4RI09IIDSEGAZzDQ15DECBxAGgI8g/0jcBZ1AJCAAM/gwAicKRQAxRHgGQz2M4UCIB0BnsFA7+JBTAIcQHEEeAaLI/C9TKEASECAZ9Bqz4FzNiaAFAR4Bq12FwojcAApCPAM+ps9shmBA0hDgGdgL/SC76UPHEACAjyTgb4GfeAAkhDgmQz2NWgjBJCEAM+k1ddkBA4gCQGeSYsROIBEBHgmrf4my6oBSJIU4Lb/1vaDtrfZvsb2irIK63at3oammEIBkKBwgNteL+kzkkYj4gxJDUmXlFVYt2v1N2gjBJAkdQqlKWnAdlNSS9KO9JKWh1ZfQ1P7CHAAxRUO8Ih4UtI/SXpc0k5JL0TEzw59n+1Ntsdsj01MTBSvtMssnMjDFAqA4lKmUIYkXSzpDZLWSRq0femh74uIzRExGhGjw8PDxSvtMgN9DU3PzmtuPnKXAqCmUqZQ3ifpNxExERH7JF0n6ffLKav7sbAxgFQpAf64pHNtt2xb0gWSxsspq/sNcElZAIlS5sDvkvR9SVslPdD+rs0l1dX1DizqQIADKKaZ8uGI+LKkL5dUy7JyYFk1plAAFMOZmJmwsDGAVAR4JovLqtFKCKAoAjyTxWXVOIgJoCgCPJPB/etiEuAAiiHAMznQRsgUCoBiCPBMaCMEkIoAz2RFky4UAGkI8Ex6erywKg9dKAAKIsAzavU1NMklZQEURIBn1OprMgIHUBgBnhELGwNIQYBnRIADSEGAZ9Tqa3IxKwCFEeAZMQIHkIIAz4gAB5CCAM+o1c8UCoDiCPCMWr2MwAEUR4BntDACn9M8K9MDKIAAz2hw8YqEnI0JoAACPKPFVXmYRgFQBAGe0eKqPBzIBFAEAZ4RI3AAKQjwjFr9jMABFEeAZ8QIHEAKAjyjxQDfO02AA+gcAZ7R4kHMqX1MoQDoHAGe0SAjcAAJkgLc9irb37f9kO1x2+8qq7DlYGDxRB7mwAEU0Ez8/D9LuiEi/sR2n6RWCTUtG4tTKHvpQgFQQOEAt328pHdL+gtJiogZSTPllLU8NHqs/mYPI3AAhaRMobxR0oSkb9q+1/bXbQ8e+ibbm2yP2R6bmJhI2Fx3GuxvMgIHUEhKgDclnSPp3yLibEl7JX3h0DdFxOaIGI2I0eHh4YTNdScWdQBQVEqAb5e0PSLuaj/+vhYCHR1o9TU0SRcKgAIKB3hEPCXpCduntZ+6QNKvSqlqGWn1MYUCoJjULpS/kXR1uwPlUUmfSC9peWn1NTiICaCQpACPiPskjZZUy7LU6mvq+cmp3GUAqCHOxMxsYQTOFAqAzhHgmQ32N7SXKRQABRDgmQ30NpkDB1AIAZ7Zwgh8VhGsTA+gMwR4ZgN9DUVI07PzuUsBUDMEeGaDixe0muZAJoDOEOCZDbCsGoCCCPDMFkfgBDiAThHgmR1Y2JgpFACdIcAzY2V6AEUR4JkN9jOFAqAYAjyzAaZQABREgGd2oI2QETiAzhDgmTECB1AUAZ4ZBzEBFEWAZ9bb6FFfo4cAB9AxArwCWv0NplAAdIwAr4BWLyvTA+gcAV4Brf4mI3AAHSPAK+CEgV49t3cmdxkAaoYAr4CRoQE9uZuFjQF0hgCvgJGhAe3c/bJm51jUAcDSEeAVMDLU0ux86OmXpnOXAqBGCPAKGBkakCRtf24ycyUA6oQAr4CRoZYkafvzzIMDWDoCvALWrVohiQAH0BkCvAL6mw2ddHy/tj/PFAqApUsOcNsN2/fa/nEZBS1XI0MtRuAAOlLGCPwKSeMlfM+yNjI0oO27GYEDWLqkALc9IunDkr5eTjnLF73gADqVOgL/qqTPS3rV1LG9yfaY7bGJiYnEzXUvesEBdKpwgNv+iKRdEbHltd4XEZsjYjQiRoeHh4turuvRCw6gUykj8PMkXWT7t5K+K+l821eVUtUyRC84gE4VDvCI+GJEjETERkmXSLolIi4trbJlhl5wAJ2iD7wi6AUH0KlmGV8SEbdJuq2M71rO6AUH0AlG4BVCLziAThDgFbLYCz43H7lLAVADBHiF7O8Ff/Hl3KUAqAECvEL294IzDw5gCQjwCjnQC848OIAjI8ArhF5wAJ0gwCuEXnAAnSDAK4ZecABLRYBXzMjQAAEOYEkI8IoZGRrQjt1T9IIDOCICvGLoBQewVAR4xdALDmCpCPCKoRccwFIR4BVDLziApSLAK4ZecABLRYBXEL3gAJaCAK+g9avoBQdwZAR4BdELDmApCPAKohccwFIQ4BVELziApSDAK+hAgNOJAuDVEeAVtG4VI3AAR0aAV9CK3oZOPI5ecACvjQCvKC4rC+BICPCK2rC6pUd27VEErYQADo8Ar6h3v3lYu16a1tbHn89dCoCKIsAr6gNvPVn9zR5df9+O3KUAqKjCAW77FNu32h63/aDtK8osbLlb2d/UBaefqJ88sFOzc/O5ywFQQSkj8FlJn4uI0yWdK+ly228ppyxI0kVnrtMze2b0i0efy10KgAoqHOARsTMitrbvvyRpXNL6sgqD9J7TTtTK/qauv//J3KUAqKBS5sBtb5R0tqS7DvPaJttjtscmJibK2NyysaK3oQ+89ST957anND07l7scABWTHOC2V0q6VtJnI+LFQ1+PiM0RMRoRo8PDw6mbW3YuOnOdXnp5Vj9/mB8/AK+UFOC2e7UQ3ldHxHXllISDnfemNRpq9er6++lGAfBKKV0olvQNSeMR8ZXySsLBehs9+tDb1urm8V2anJnNXQ6ACkkZgZ8n6WOSzrd9X/vPh0qqCwe56Mx1mto3p5t+9XTuUgBUSLPoByPiDkkusRa8irdvXK2Tj1+hH92/QxefRaMPgAWciVkDPT3WR353rX7+PxPaPTmTuxwAFUGA18RFZ63TvrnQjQ8+lbsUABVBgNfE29afoI2va9GNAmA/ArwmbOvis9brv//3Wf30gZ25ywFQAQR4jWx69xt1zoYhfeaae3XLQ3SkAMsdAV4jg/1NffMTb9fpa4/XX121VXf8+pncJQHIiACvmeNX9Oo7l71Db1wzqE995x7d9eizuUsCkAkBXkNDg33690++U+tXDeiyb92je1m1B1iWCPCaGj6uX1d/6ly9bmW//vzKu/WVnz2sx57dm7ssAMcQAV5jJ5+wQv/xl+/U2RuG9C+3PqI/+Mfb9Gdfu1Pfu+cJ7ZnmuilAtyt8Kj2qYWSope9c9g7tfGFK1219Utdu2a7PX/tL/d312/Tmk47ThtUtbVjd0utf19Ipq1tas7Jfrb6GVvY31eprqq/JbzhQVwR4l1h7woAuf++b9On3/I7ufWK3fnT/Dj2ya48eePIF3bDtKc3Ox2E/19foUX+zR82G1ejpUW/DavRYzR7LtuyFC9707L+/cLvI9mEviGOukgO8wtcu/T2dsrpV6ncS4F3Gts7ZMKRzNgztf252bl47X3hZjz07qd1TM5qcntOe6VlNzsxqz/ScpmfnNDcf2jcXmpuf1+xcaHY+FJIiQhFSKDQ/v3Ab7d+ChdcPV8XhfyyA5azZKH9UQ4AvA81Gj05Z3Sr91x9AXkyAAkBNEeAAUFMEOADUFAEOADVFgANATRHgAFBTBDgA1BQBDgA1RYADQE0R4ABQUwQ4ANQUAQ4ANUWAA0BNJQW47QttP2z7EdtfKKsoAMCRFQ5w2w1J/yrpg5LeIumjtt9SVmEAgNeWcj3wd0h6JCIelSTb35V0saRflVHYwb70gwd092+eK/trs3i1lXEAoFMpAb5e0hMHPd4u6Z2Hvsn2JkmbJGnDhg2FNrRu1YBOPWlloc9WzZtPPk4fPGNt7jIAdIGUAD/c+kD/b3gZEZslbZak0dHRQsPPy9/7piIfA4CulnIQc7ukUw56PCJpR1o5AIClSgnweySdavsNtvskXSLp+nLKAgAcSeEplIiYtf3Xkm6U1JB0ZUQ8WFplAIDXlLQqfUT8VNJPS6oFANABzsQEgJoiwAGgpghwAKgpAhwAasoRx+7UbtsTkh4r+PE1kp4psZycumVfumU/JPalqrplX1L34/URMXzok8c0wFPYHouI0dx1lKFb9qVb9kNiX6qqW/blaO0HUygAUFMEOADUVJ0CfHPuAkrULfvSLfshsS9V1S37clT2ozZz4ACAV6rTCBwAcBACHABqqlYBbvtPbT9oe9527VqLumURaNtX2t5le1vuWlLZPsX2rbbH23+3rshdU1G2V9i+2/b97X35+9w1pbDdsH2v7R/nriWF7d/afsD2fbbHyvzuWgW4pG2S/ljS7bkL6VSXLQL9LUkX5i6iJLOSPhcRp0s6V9LlNf73Mi3p/Ig4U9JZki60fW7mmlJcIWk8dxEleW9EnFV2L3itAjwixiPi4dx1FLR/EeiImJG0uAh07UTE7ZK6YpXpiNgZEVvb91/SQmCsz1tVMbFgT/thb/tPLbsUbI9I+rCkr+eupcpqFeA1d7hFoGsZFN3K9kZJZ0u6K28lxbWnHe6TtEvSTRFR1335qqTPS5rPXUgJQtLPbG9pL/JemqQFHY4G2/8l6eTDvPSliPjhsa6nREtaBBp52F4p6VpJn42IF3PXU1REzEk6y/YqST+wfUZE1OpYhe2PSNoVEVtsvyd3PSU4LyJ22D5R0k22H2r/X2yyygV4RLwvdw1HCYtAV5TtXi2E99URcV3uesoQEbtt36aFYxW1CnBJ50m6yPaHJK2QdLztqyLi0sx1FRIRO9q3u2z/QAvTqaUEOFMoxw6LQFeQbUv6hqTxiPhK7npS2B5uj7xle0DS+yQ9lLeqzkXEFyNiJCI2auG/k1vqGt62B20ft3hf0gdU4g9qrQLc9h/Z3i7pXZJ+YvvG3DUtVUTMSlpcBHpc0vfqugi07Wsk3SnpNNvbbX8yd00JzpP0MUnnt9u87muP/OporaRbbf9SCwOGmyKi1i14XeAkSXfYvl/S3ZJ+EhE3lPXlnEoPADVVqxE4AOAAAhwAaooAB4CaIsABoKYIcACoKQIcAGqKAAeAmvo/OFRRuSlt5pwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot_step_response(Gc)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.4" } }, "nbformat": 4, "nbformat_minor": 4 }